SIZE-STRUCTURED POPULATIONS: IMMIGRATION, 
(BI) STABILITY AND THE NET GROWTH RATE 



JOZSEF Z. FARKAS 



Abstract. We consider a class of physiologically structured population 
models, a first order nonlinear partial differential equation equipped with 
a nonlocal boundary condition, with a constant external inflow of indi- 
viduals. We prove that the linearised system is governed by a quasicon- 
traction semigroup. We also establish that linear stability of equilibrium 
solutions is governed by a generalized net reproduction function. In a 
special case of the model ingredients we discuss the nonlinear dynamics 
of the system when the spectral bound of the linearised operator equals 
zero, i.e. when linearisation does not decide stability. This allows us 
to demonstrate, through a concrete example, how immigration might be 
beneficial to the population. In particular, we show that from a nonlin- 
early unstable positive equilibrium a linearly stable and unstable pair of 
equilibria bifurcates. In fact, the linearised system exhibits bistability, 
for a certain range of values of the external inflow, induced potentially by 
Allee-effect. 

1. Introduction 

Interest in the modelling and understanding of the dynamics of biological 
populations is old. Classical, ordinary differential equation models assume ho- 
mogeneity of individuals within population classes, and therefore involve equa- 
tions for total population sizes. However, individuals in biological populations 
differ in their physiological characteristics. Therefore vital rates, such as those 
of birth, death and development, vary amongst individuals. Structured popu- 
lation modelling tries to account for these differences. It has played a central 
role in mathematical population dynamics in the past decades, being applied 
to amongst other systems, fish, water flea, and cell populations. The mathe- 
matical modelling of tumour growth and cancer treatment, human demography 
and the epidemiology of infectious diseases may also involve structured pop- 
ulations. The first step when modelling a population may be to consider a 
variable that allows one to divide the population into internally homogeneous 
sub-populations in order to describe the dynamics as the interaction of these 
groups. This structuring variable may be discrete [8] or continuous [20l [2T1 [26] 
and may be based on age, size, stage of development or life cycle. Tradition- 
ally, (continuous) structured population models have been formulated as par- 
tial differential equations, starting with the pioneering (first nonlinear) model 
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of Gurtin & MacCamy pjj]. Following the lead of [22J, [23] [26] we successfully 
applied the spectral theory of compact /positive operators to formulate biologi- 
cally interpretable conditions for the linear stability/instability of equilibria for 
age and size-structured population models (see [T2] EH HH [15l [16] El HB] ) ■ Our 
method allowed us to formulate very general and elegant stability/instability 
conditions for scramble competition models in terms of certain associated re- 
production functions. We shall note that previously, in 24 , Rorres formulated 
linear stability/instability conditions in terms of the derivative of a net repro- 
duction function for positive equilibrium solutions of an age-structured model 
in the special case when only the fertility depends on the standing population 
size. We shall also point out that Rorres' analysis in [24] was not rigorous. 

Although traditionally structured population models have been formulated 
as partial differential equations for population densities, the recent unified ap- 
proach of Diekmann et al. , making use of the rich theory of delay and integral 
equations, has been resulted in significant results. Most notably the Principle 
of Linearized Stability has been proven in [9] QJ5] for a wide class of physiologi- 
cally structured population models formulated as delay equations (or abstract 
integral equations). 

In the recent papers [T2] [TB] we investigated analytically the effects of the 
introduction of a positive influx of newborns on the linearised dynamical be- 
haviour of age and size-structured scramble competition models. In our main 
result we were able to relax our linear stability conditions given previously in 
12, 14 for models without inflow. Motivated by our new stability results we 
introduced the concept of a net growth rate and formulated a very general con- 
jecture in [13] , One of our goals in the present work is to show that our stability 
results in [14] can be extended analogously in full generality to size-structured 
scramble competition models with (constant) inflow at birth using the recently 
introduced concept of the net growth rate. 

In particular, in the present paper we consider the following nonlinear first 
order partial differential equation. 
(1) 

Pt (s,t)+(j(s,P(t))p(s,t)) s = -iJ,(s,P(t))p{s,t), t>0,se(0,m), m<oo. 

Equation (UJ) is equipped with the following nonlinear nonlocal boundary con- 
dition representing the influx of individuals into the population 

(2) 7 (0,P(t))p(0,t) = C+ / (3(s,P(t))p(s,t)ds, t>0, 0<C<oo, 

Jo 

and is subject to the initial condition 

(3) p(s,0)=:p (s), se(Q,m). 

Here, the function p = p(s, t) denotes the density of individuals of size s at 
time t, hence P(t) = J Q p(s, i) ds gives the total population quantity at time t. 
(3, fi and 7 denote the fertility, mortality and development rates of individuals, 
respectively, and all of these vital rates depend on both the structuring vari- 
able s and on the total population size P(t). We make the following general 
assumptions on these vital rate functions: 

H, 13 e C l {[Q,m] x [0,oo)), /3,/x>0, 7 € C 2 ([0, m] x [0, 00)), 7 > 0. 
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From the biological point of view it is natural to assume a finite maximal size 
m which in turn helps to avoid mathematical difficulties regarding the spectral 
analysis of the linearised system (see e.g. |170. C represents an inflow of min- 
imal size individuals, i.e. an inflow of individuals from an external source at 
birth. There are several biological motivations for having a population inflow of 
individuals from an external source. In a discrete model the inflow might rep- 
resent artificial stocking of individuals at certain life cycle stages for example in 
an insect population. In [25 Smith treated a nonlinear stage-structured matrix 
model with stocking and he applied his results to the celebrated "LPA model 
for the flour beetle Tribolium" to discover some interesting phenomenon due to 
the stocking. In the continuous setting, in which we are interested, a natural 
example for such an inflow is the case of migratory fish populations (such as a 
number of salmon species) that lie eggs and then move on to a different habi- 
tat. Later on the newly hatched fish join a different local population. Another 
example underscoring the relevance of this model is the case of fisheries where 
stocking of newborn (minimal size) fish is practiced. In this scenario we can 
regard the inflow C as a control parameter in our system. 

Model (p}-© was previously studied in [4] by Calsina and Saldana. They es- 
tablished global well-posedness and also investigated the asymptotic behaviour 
of solutions. In particular, they showed that under certain conditions on the 
vital rates H])-© is governed by a strongly continuous semigroup which ad- 
mits a compact global attractor. Most notably, they obtained conditions on 
the model ingredients which guarantee convergence of solutions to a stationary 
distribution when the total population tends to a constant. In contrast, our 
paper focuses on the local asymptotic stability of equilibrium solutions of sys- 
tem ©-(GO) with particular regards to the effects of the population inflow on 
the dynamical behaviour of the system. First, we will prove the existence of 
a quasicontraction semigroup describing the evolution of solutions of the lin- 
earised system. We also discuss regularity properties of the semigroup which 
are crucial for our stability studies. Next we will give stability/instability crite- 
rions for positive equilibrium solutions of (jT|)-([3]) . Then, utilizing the stability 
results and a "naive" approach, we will investigate the nonlinear dynamics of 
the model in a special case of model ingredients. In particular we will treat the 
case when, in the abscence of inflow, the strictly dominant eigenvalue of the 
linearised operator is zero. We will then show, through a concrete example, 
that the introduction of the inflow induces bistability in the system. Our claims 
will be supported by numerical examples, as well. 



2. Equilibrium solutions, semigroup existence and linear stability 

In this section we investigate the existence and linear stability of equilibrium 
solutions of (JTJ)— (JSJ) - We note that (p}-© admits the trivial solution if and only 
if C = 0. When we solve H])-© for a positive stationary solution p„ we arrive 
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at the following equations: 
(4) 

p*(s) = p*(0)7r(s, P*), 
(5) 

pm pm 

-f(0,P,)p,{0) = C+ (3(s,P*)p4s)ds = C+p*{0) (3(s,P*)ir(s,P*)ds, 
Jo Jo 

where 

(6) 7 r( S ,P,) = exp|-y o drj. 

Thus, for given vital rates /?, /u, 7 and inflow C the function is a positive 
stationary solution of dU)- (J3j) if it is determined by 

(7) p.(«) = r^'iV. , «e(0,m), 

Jo As-,P*)as 

with the positive total population quantity 

/■771 

(8) P* = / p*(s)ds, 

Jo 

which satisfies the equation 

(9) 7(0,P*)P*=C/ ir(s,P*)ds + P* f3{ S) P*)ir(s,P*)ds. 

Jo Jo 

Then, if for fixed C G [0, 00) we define the net growth rate Qc f° r P £ (0, 00) 
as follows: 

(10) Q C (P) = -TjrprfC P^ 1 / 7r(s,P)ds+ P(s,P)Tr(s,P)ds), 

7(0, P)V 7 Jo 7 

it is easily shown that positive equilibrium solutions of J]}-© are in a one-to- 
one correspondence with positive solutions P* of equation 

(11) Qc(P) = i. 

Next we note that 

1 r^ PMSjP)ds= fM^Lf^L 



7 (o,p)7 , v , , y 7 ( s ,p) ^ y 7 ( r ,p) 



(12) 



y /3(a, P) exp j - y /i(r, P) dr j da, 



by changing variables from size s to age a using the relation 

(13) £=7(«(a),.)- 
Hence the function P, defined by 

(14) i2(P) = _J_/" /3( s ,P)7r( S ,P)d S = g (P), 

71 u j Jo 

is the well-known inherent net reproduction function, i.e. the expected number 
of newborns to be produced by an individual in her lifetime when the population 
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size is P. We also notice that CP 1 is the per capita inflow and by (fT3|) wc 
have 



(15) 



where L is in turn the expected lifetime of an individual. Therefore (|10[) can 
be rewritten as 

(16) Q C (P) = J^ (CP' 1 + /3(a,P))exp j- J n(r,P)dr\ da. 

It is straightforward to verify that the rather natural assumption 

(17) lim P(P) = 0, 

P^oo 

yields a sufficient condition for the existence of at least one positive solution of 
equation (|11[) and in turn for the existence of at least one positive equilibrium 
solution of (P)-© via formula (J7]). 

Given a positive stationary solution of system ([l|)-((3]), we introduce the 
perturbation u = u(a,t) of p by making the ansatz p = u + p*. Then we 
are using Taylor series expansions of the vital rates to arrive at the linearised 
problem 

u t (s,t) +j(s,P*)u s (s,t) + (7 s (s,P») + /j,(s,P*)) u(s,t) 
(18) 

+ (7 s p(s, P*)p»(s) + fJ.p{s, P*)p*(s) + 7p(s, P*)p*'(s)) U(t) = 0, 

(19) 

ln+ , r ( ^ P *) - 7p(0, P.)P*(0) + C M*", P)p*M dr \ 

» ( o,t) = y o ^ ^-^ ju(.,t)* 

where we have set 

(20) t7(t) = / u(s,t)ds. 

Jo 

Eqs. (|18 p - (fT9|) are accompanied by the initial condition 

(21) u(s,0) =u (s). 

Linear semigroup theory proved to be a usefool tool to study asymptotic be- 
haviour of solutions of similar types of physiologically structured population 
models, see e.g. [HI HH [15j [16j H3 HH] ■ The first crucial step in such an analy- 
sis is to establish the existence of the governing linear semigroup. In [14| . IT5 | [16] 
we used the Desch-Schappacher perturbation theorem to establish the existence 
of the governing semigroup. More recently, in [171 118j we utilised dissipativity 
calculations to establish the existence of a governing quasicontraction semi- 
group. Here we present another approach, which is a little bit of a mix of the 
above two but uses the Trotter product formula, and it is potentially more 
useful since further regularity properties of the semigroup, such as positivity, 
may be immediately established. 
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To this end, we cast the linearised system (jl8p -([2l j) in the form of an abstract 
Cauchy problem on the biologically relevant state space X = L 1 (0,to) as fol- 
lows: 



(22) 
where 



^u = (A + B + C)u, u(0) 



Au = (7(-,P*)u) with domain Dom(^) = {u 6 W 1A (0, m) I u(0) = <f>(i 

as 1 

Sit = —//(•, P*) u on ,Y, 

Cu = - (7 s p(-,P*)p» + /Up(-,P*)f>* + 7p(',P*)pQ / u(s)ds 



$(u) 



tt(s) as on A 7 , 

' ) - 7P (0, -P* )P« (0) + Jo"' /?p (r, ) P* (r) dr ' 

7(0, P*) 



m(s) ds on A" 



Our objective here to establish that A + B + C is a generator of a positive 
semigroup is to apply the following result (see e.g. Sect. 5 in [11]). which is a 
variant of the Trotter product formula. 

Theorem 1. Let T> and J- be generators of the strongly continuous semi- 
groups {T(t)} t>0 and {S(t)} t>0> respectively, satisfying 



(23) 



T - 5 - 



< M exp{wi} for all (>0,neN, 



for some constants M > 1, u> £ R. Assume that y = Dom(T>) n Dom(J-) is 
dense in X and for some \q > uj (Ao — T> — T)y is also dense in X. Then the 
closure ofD + J- generates a strongly continuous semigroup {U(t)} t >o given by 
the formula 



(24) 



U{t)x = lim 



T - S - 



x € X. 



Our first goal is to establish a condition which guarantees that A generates 
a positive quasicontraction semigroup. We recall (see e.g. [3]) that the linear 
operator C is called dispersive if for every y € Dom(£) C y there exists a 
4> y S y* with < (j>y, \\(j> y \\y < 1 and (y, 4> y )- = \ \y + \\y such that 

(A/,<^>- <o, 

where (•,•)- stands for the natural pairing between elements of y and its dual 
y*. Our aim is to apply a variant of the Lumer-Phillips Theorem, see e.g. 

mm- 

Theorem 2. The operator A generates a strongly continuous and positive 
quasicontraction semigroup {TA(t)}t>o of bounded linear operators on X if the 
following condition holds true: 

(25) /3(s,P*)- 7P (0,P*)p*(0)+ / l3 P (r,P*)p*(r)dr > 0, for a€[0,m]. 
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Proof. We consider the operator A — uiX for some uj £ ML For a fixed u £ X 
we define (j> u £ A"* by <^«(s) = 1 if u(s) > and <j> u {s) = if u(s) < 0. For 
u £ Dom(^l) we obtain the following estimate: 



((A - lu!)u, (f) u )- 



T (7(s, P*)u(s)) Xu+ (s) ds-uj 
as 



u{s)Xu+(s)ds 



< - 7 (m, P*)«(m) X „+ (m) + 7 (0, P*)u(0)x„+ (0) - u\ \u + \ | a 



(26) 



= &||w H 
<0, 



■«||u+||i 



/3p(?-, P*)p*(r) dr > < oo 



for some uj > b large enough, where 

6< sup j/3(s,P*)-7 F (0,P*)p*(0) + 

s£[0,m] I 

Hence we have showed that for some u> £ K large enough the operator — cjI 
is dispersive. 

To obtain the range condition, i.e. that Rg(AZ — A) = X we obtain the 
unique solution of the equation 

(27) (XI - A)u = h, 

for A large enough and h £ L 1 (0, m) as follows: 

A + 7s (y,P*) 



it(s) = exp 



(28) x ( u(0) + 

Then we have 



exp 



dy 



v A + 7s (r,P») 



\u'(s)\ < 



h(s) 



(29) 



7(s,P 
exp 



|«(0)| 



7(r, P, 
A + 7s(s,P*) 



d.r 



■ dy 



A + 7s (r,P») 



exp 



j(y,P*) 

A + 7s (y,P,) 



7(1/, -P. 



dy 



dr 



i(y,P*) 



dy. 



That is u £ W 1 ' 1 (0,m) and the range condition is satisfied and the proof is 
now completed on the grounds of Theorem 1.2 in Section C-II in [3]. □ 



Theorem 3. Assume that condition (|25|) and the following condition hold 
true 

(30) P*(s)>0, for s£[0,m}. 

Then, A + B + C generates a positive quasicontr action semigroup. 

Proof. If condition ([3H| holds true then C is a positive operator. Since B + C 
is bounded and linear and B is a multiplication operator it follows that B + C is 
a generator of a positive semigroup if condition (|30[) holds true. Finally, since 
both A and B + C are generators of positive semigroups, Theorem 1 shows that 
A + B + C is a generator of a positive quasicontraction semigroup if conditions 
13 and ((201) bold true. □ 
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Next we would like to establish conditions which guarantee that the equilib- 
rium solution is linearly asymptotically (un)stable. In particular, we would 
like to give a complete generalization of our stability results in [14] for model 
(H}-© with boundary immigration. To formulate stability/instability theo- 
rems one usually deduces a characteristic function (if possible!) which roots 
are the eigenvalues of the linearised operator A + B + C. Then, for instability 
one only needs to establish the existence of a positive root of this characteristic 
function. For stability however, one also needs to guarantee that the growth 
bound of the semigroup, and hence stability, is governed by a real eigenvalue. 

Theorem 4. The spectrum of A + B + C can contain only isolated eigen- 
values of finite algebraic multiplicity. 

Proof. Our aim is to show that the resolvent operator R(X, A + B + C) is 
compact. Since B + C is bounded, it is enough to show that R(X, A) is compact. 
We obtain the solution of the resolvent equation 

(31) (X-A)- 1 f = u 

as 







We note that, for A large enough the solution u of the resolvent equation 
(f3"Tj) belongs to M /1,1 (0,m). That is for A large enough the resolvent operator 
{XT — A)^ 1 is a bounded linear mapping from L l (0, m) into W ' (0, m). Since 
Fy 1,1 (0, m) is compactly embedded in £ 1 (0, m) by the Rellich-Kondrachov the- 
orem (see e.g. [2]), the claim follows. □ 

We are now able to formulate our main stability result. 

Theorem 5. Let be a non-trivial stationary solution with corresponding 
population quantity . If 

(33) Q'c(P*) > 

then is linearly unstable. (Note that Q' c stands for the derivative of the 
function Qc with respect to P.) On the other hand, suppose that conditions 
(f2"S")) - (|3TH) and the following condition hold true: 



P*(s) 
o Jo 



j s p(r,P*) + fJ.p(r,P*) 
l(r,P*) 
Mr,P*){ls(r,P*)+Kr,P*)) 



(34) „ 

T[r,P*) 

Then, p* is linearly asymptotically stable if 

(35) Q' C {P*) < 0- 



dr ds > — 1 . 
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Proof. It is easily shown (see e.g. [14]) that the eigenvalue equation 

(36) {A + B + C - XI) u = 0, u(0) = $(u), 

admits a solution A G C with a corresponding non-zero eigenvector u, if and 
only if A satisfies the following characteristic equation: 

i = K(X) = J /(A, s )b( s ) ds (l - £ l /(A, s) J' r x (A, y)^fry d y ds 

+ / m /(A, S ) f S f-\X,y)-EM-dyd s 

Jo Jo j{y,P*) 

p rn p rn ps / \ 

(37) +/ f(X,s)ds f(X,s)b(s) f-i(\,y)-^±-dyds, 

Jo Jo Jo 7{y,P*) 

where we defined the following functions: 

/(A,s) =cxp|-y — - dyj, se[0,m],AeC, 

7(0, P*) 

It follows from the growth behaviour of the function / introduced above that 
the characteristic function satisfies 

lim K(X) = 0. 

A— > + OC 

Straightforward claculations yield: 

p m p rn ps r \ 

K(0) = - TT(s,P,)b{s)ds 7r(«,P.) / ^—TTT^r-^dyd 3 



o n{y,P*h{y,P*) 

pm pm ps f \ 

+ / TT(s,P*)b( S )d S + / tt( S ,P*) / —J^tL—dyds 

Jo Jo Jo ^{y,P*n{y,P*) 

pm pm ps t \ 

+ TT(s,P*)ds n(s,P*)b(s) . ?\ [y ^. dyds 



o 7r (y>P*)7(2/,P*) 
ir(s,P*)b(s)ds + P„ / TT P (s,P*)b(s)ds+p*(0) / Tr P (s,P*)ds 



(38) -p*(0) / ir P (s,P*)ds / 7r(s,P»)&(s)ds, 
Jo Jo 

where 7rp stand for the derivative of the function tt with respect to P and we 
used the relation 



P* =p»(0) / Tr(s,P*)ds. 
Jo 

On the other hand, we have 

^ (/ o m 7T( S ,P)d S ) \ P=P , X C 7 (0,P.)P* 



Q^(P*) = P'(P») 
(39) 



( 7 (0,P,)P*) 2 
C( 7P (0, P,)P, + 7 (0, P*)) /*(«, P.) da 



( 7 (0, P*)P*) 2 
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and 

fl'(P*) 7 (0,P*) 

r rn 

(40) =/ (p P (a,P m )ir(s,P.) + P(8,P*)np(s,P.))da-R(P*)'Yp(0,P.). 
Jo 

Further algebraic manipulations of equations (f39|) -(j40 |) yield: 
P*Q' C (P*) = ^JT ) Pp(s, P*Ms, P*)ds + jf P*)7r F (s, P.) ds^j 

7(0, P.) 7(0,P»)7 

1 ^ Jo ' V7 2 (0,P*) P*7(0,P*)7 

Comparing equations (|38p and (|4ip we arrive at 
(42) X(0) - P*Q' C (P») + 1. 

Therefore, if condition (|3"3")l holds true then from the Intermediate Value Theo- 
rem it follows that there exists a positive eigenvalue and the stationary solution 
is linearly unstable. 

Next we observe that condition (|34|) guarantees that the characteristic func- 
tion K{X) defined in (|37p is monotone decreasing for A > 0. Therefore, equa- 
tion (j42|) shows that A + B + C does not have a non-negative real eigenvalue if 
Q' C (P*) < holds. Next we recall that for the strongly continuous semigroup 
{T(t)} t>0 with generator A + B + C we have 

ujq = max {s(A + B + C),uj ess (A + B + C)} , 

see e.g. [TT], where uj ess stands for the essential growth bound of the semigroup. 
However, for a positive semigroup the spectral bound belongs to the spectrum 
unless the spectrum is empty. Theorem 4 implies that the essential spectrum 
of the generator A + B + C is empty, therefore the growth bound is determined 
by a leading eigenvalue unless the spectrum is empty. In either case, we obtain 
for the growth bound loq of the linearised semigroup 

uj = s(A + B + C) < 0, 

and the proof of the theorem is completed. □ 



3. Nonlinear dynamics and bistability 

In [14j our aim was to formulate linear stability/instability conditions for 
the model ([I])-© in the abscence of inflow. Therefore, the case i?'(P*) = 
Qo(P») = was left completely untreated. In this case, zero is an eigenvalue 
of the linearised operator, see equation (|42|) . If the positivity conditions (f25|) - 
(|30p also hold true, then zero is in fact a strictly dominant eigenvalue, hence 
linearisation does not decide stability. The goal of the present section is to 
consider this interesting and rather difficult case. To our knowledge there are 
a very few results in the literature (see [7]) which treat the case when the 
spectral bound of the linearised operator equals zero. We point out that our 
discussion will be less rigorous and some of the details need to be elaborated. 
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For basic concepts and results used throughout this chapter we refer the reader 

to |a ED. 

For the rest of the section we assume that conditions (p25")) , (13TJ)) and (|3"4")) 
hold true. The semigroup {T(t)}t>o generated by the operator A + B + C 
is eventually compact. The spectrum of the linearised operator A + B + C, 
consists of isolated eigenvalues of finite (algebraic) multiplicity. Therefore we 
may split the spectrum a {A + B + C) into the union of two disjoint closed 
subsets as follows: 

(43) a{A + B + C) = {0}U (cr(A + B + C) \ {0} 

Moreover, X splits into the direct sum of two {T(i) } t >o-invariant closed sub- 
spaces: y and 2. In fact, there exists a spectral decomposition of X with 
respect to (|4"3"f . This is because the set {0} in is trivially bounded, see 
e.g. [3]. Note that, in general, the set a (A + B + C) \ { 0} is non-empty, hence 
the spectral decomposition is not trivial, but the formulation of precise con- 
ditions is left for future work. Moreover, we may determine the subspace y 
explicitly. Straightforward calculations show that the eigenvalue equation 

(44) (A + B + C)U = 

subject to 
(45) 

TTfn , f m ( g(g, P*) - 7p(0, P*)P*(0) + Jo" Mr, P*)p*(r) dr \ TU , , 

m= l { w^) ) u{s)ds > 

has a one-dimensional solution space spanned by the function 

/ 1 + C F(s) f n s §4 dr ds r G(r) \ 

(46) F(s)[ J ° ; w\ d A e^(0,m), 

I Jo F{s)ds Jo F(r) J 

where 

(47) F{s)=ev{- n-W+^'V *), se[0,m], 

I Jo l(r,P*) J 

and 

(48) G(r) =-y aP (r,P,)p*(r) + fjL P (r, P,)p,(r) + 7P (r, P*K(r), r e [0,m]. 
Hence the one-dimensional subspace y of X given by 

11 \ { '{ rm<h J« F ir) 

which is (A + B + C)-invariant, is the centre eigenspace. We note that the set 
a(A + B + C)\ {0} is contained in the half plane Hg ={zeC : Re(z) < 6} for 
some 6 < 0. This, and the existence of the spectral decomposition implies that 
the semigroup {S(t)}t>o, which is the restriction of {T(i)} t >o to Z is uniformly 
exponentially stable (see e.g. [3]). There exists a centre manifold (see e.g. [6]), 
but it cannot be determined explicitly since the eigenvalues of A + B + C (and 
in turn the corresponding eigenspaces) are not explicitly available (except the 
zero eigenvalue). 
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In the rest of the section we restrict our attention to the special case of vital 
rates: 



R'{P*)=j^ Ms,P*)expj-^ 



dr > ds = 0, 



(50) M = /*(*). 7 = 7(«), P = P{s,P), s£ [0,77i], Pe[0,oo), 
and we further assume that 

7 a (r) + M (r) 
7(r) 

(51) #'(P.) = jT P PP (s, P m ) exp |- jf 7 " ( ^y W drj + 0. 

In this scenario, equation (fTJ) is linear, and the only nonlinearity arises in the 
boundary condition ([2]). Since linearisation does not decide stability we take 
into account the quadratic term when using Taylor series expansion for the 
fertility function (3. We obtain the nonlinear problem: 

jU = (A + B)U, 
Dom(^) = 

U€W 1 <H0,m)\U(O)= r^l^lu(s)ds + U 2 C M "' P * ] W *) ds 
Jo 7(0) Jo 7(0) 

(52) 

p rn 

U= U{s)ds. 



When we look to obtain solutions of the equation 

(53) (A + B)U = s, eeR, 
we arrive at 

(54) = e/ /?(s,P*)F(s) / -—drds + U 2 f3 PP (s, P*)p*(s) ds, 

Jo Jo F { r ) Jo 

an equation for U. Existence of a (real) solution of (|54|) is shown to be nec- 
cessary and sufficient for the existence of a non-trivial solution of (|53ll . This 
shows that (|53|) admits solutions only for e > or for e < depending on 
the sign of R"(P*). This means that the flow generated by the linear operator 
A + B is monotone in a neighbourhood of U = 0, which in turn implies that 
the equilibrium is unstable. 

This claim will be confirmed numerically in case of an example. When we 
make the following choice of vital rates: 

P 2 e- p se~ s + 0.5P 2 e- p 
3e" 2 - 2e- 8 - 13e" 
system ([T])-([3]) admits for C = a unique positive equilibrium solution 



( 55 ) M = l, 7S1, /? = o-^,,-,,^ ^[0,6] 



(56) Pa; ( s ) = 2-^— -, se [0,6], 

1 — e D 

with total population size P* = 2. It is straightforward to confirm that con- 
ditions ([23, (001) and ([31]) hold true and Q (P») = P'(P*) = 0, Q '(P*) = 
R"(P*) < since /3 PP (s,P») < holds for s e [0,6]. For C* = we have 
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Qo(0) = < 1 which implies that the trivial equilibrium is locally asymp- 
totically stable (see e.g. [33]). The positive equilibrium p* given by (l56|) is 
nonlinearly unstable. 

When we introduce the inflow, at C = a pair of equilibria, denoted by pi 
and pi , bifurcates from the positive equilibrium p» . We may assume that 



(57) 



P*= / pl(8)d8<P*= / pl(s)ds. 

Jo Jo 



Then, it is easily shown that Q' C (P%) < holds, and conditions (|25j) . p0|) 
and (|34|) are satisfied. Hence the equilibrium pi is linearly stable while pi, for 
which Q'c(P-t) > holds, is linearly unstable by TheoremEl The third positive 
equilibrium, which is "connected" to the trivial one, satisfies also the stability 
conditions of Theorem hence it remains locally asymptotically stable. The 
system is bistable for C G (0, C») for some C* > . This bistability is induced 
potentially by an Allee-type effect. When C reaches the critical value C* the 
system undergoes another bifurcation. This is illustrated in Figure 1. 




FIGURE 1 . Schematic bifurcation diagram of equilibrium pop- 
ulation sizes 



In Figure 2. we present the results of some numerical simulations. We 
used the hier-community program which solves equations ( l)-(3) using a finite 
difference scheme approximation, see [T]. It is developed by A. S. Ackleh and 
his collaborators and it is available at Ackleh's web page. 
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Figure 2. In i/ie /irsf column we plotted the evolution of 
the solution and its integral starting with the initial condition 
Po(x) — ^ 9 ^_ti , and the graph of the net reproduction rate R, 
respectively. (Here x replaces the structuring variable s.) In 
the second column we plotted again the evolution of the solution 
and its integral starting with the initial population distribution 
Po(x) = ° 1 7 ^ e _ 6 , and the graph of net growth rate Q, in the 
case of C = 0.2. 

4. Concluding remarks 

In the present manuscript we employed semigroup and spectral methods to 
investigate the linearised dynamical behaviour of a size-structured population 
model. The model we treated here is an extension of the one analysed pre- 
viously in [14J . In particular our goal was to study the effects of migration 
of individuals from an external source. This question is biologically relevant. 
Migration was modelled through a constant inflow of minimal size individuals. 
The assumption of migration at only one state might seem unsatisfactory from 
the biological point of view, nevertheless we managed to give a complete gen- 
eralization of our stability/instability results established previously in [14j for 
a model without migration. In fact, in our main result we obtained biologically 
relevant stability/instability criteria by defining an appropriate net reproduc- 
tion function of the standing population. In my opinion this result is quite 
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interesting and it shall have an impact when studying other physiologically 
structured scramble competition models. 

In the second half of the present work, we were mainly concerned with the 
weakly nonlinear dynamics of the model. The bifurcation diagram in Figure 
1 illustrates the dynamic behaviour of the system in case of migration. We 
also gave an example when the introduction of the inflow stabilizes a weakly 
unstable steady state. Although in this part we treated a quite special case of 
model ingredients the analysis may well be applicable to more general scenarios. 
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